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The formation of correlations due to collisions in an interacting nucleonic system 
is investigated shortly after a disturbance. Results from one-time kinetic equations 
are compared with the Kadanoff and Baym two-time equation with collisions in- 
cluded in second order Born approximation. A reasonable agreement is found for 
a proposed approximation of the memory effects by a finite duration of collisions. 
The formation of correlations and the build up time is calculated analytically for 
the high temperature and the low temperature limit. This translates into a time 
dependent increase of the effective temperature on time scales which interfere with 
standard fire ball scenarios of heavy ion collisions. The consequences of the for- 
mation of correlations on the two- particle interferometry are investigated and it 
is found that standard extracted lifetimes should be corrected downwards. 



1 Introduction 

To extract information about the space- time character of particle- emitting 
sources is one of currently much investigated tasks an. The size of stars are 
found to be measurable by coincident^ photons d. In laboratory this method 
is used to analyse heavy ion reactions □. The two-particle correlation function 
carries the information about nuclear interaction and correlations in the emit- 
ting source Ij'LI and is given by the ratio of two - and single particle emission 
probability 



R{P = Pi+P2,q= {pi - P2)/2) + 1 - 



n(pi,P2) 
n(pi)n(p2) 



(1) 



which can alternatively be expressed by cross sectionsQ. This two-particle cor- 
relation function can be written in terms of the the two-particle wave function 
ip{q,r) 



R{P,q) + \ = J d^rFpir)\ipiq,r)\' 



(2) 



whiJe the average of the single particle Wigner functions of the source is given 

by&y 

p ^ I d'R /(^/2, R + r/2, t>) /(P/2, R - r/2, t>) 
{Jd3Rf{P/2,R,t>)f 

t> 

Here f{p, r,t >) = J dtg{p, r — p{t^ — t)/ni, t) is the phase-space Wigner dis- 
tribution of particles of momenta p at position r at some time after emission 
with probability g. A rigorous derivatiori^of this formula with the discussion of 
the necessary neglects has been given inid for the case of initially uncorrelated 
particles r)p oc 1 ±cos{rq) corresponding to bosons/fermions. Eq. (||) is 
well suited for simulation of heavy ion reactions where the one-time distribution 
function /(r,p, i) is then determined by solving appropriate kinetic equations 
(BUU) or simulating equation of motions (QMD). While most current codes 
rely on the quasiclassical Boltzmann equation including Pauliblocking effects, 
a quantum two-time theory for the time-evolution of real time Green's func- 
tions G(r,p, i,t') has been developed using the Schwinger-Keldysh formalism 
already 30 years ago. The quantum image of the classical Boltamann equation 
is usually referred to as the Kadanoff-Baym (KB) equations 113. These equa- 
tions have often been considered too complicated to sohie numerically in the 
past. However, several numerical applications exist nowtil. 

These kinetic equations describe different relaxation stages. During the 
very fast first staffp.. correlations imposed by the initial preparation of the sys- 
tem are decaying These are contained in off-shell or dephasing processes 
described by two-time propagators. During this stage of relaxation the quasi- 
particle picture is establishedliil'ta. After this very fast process the second state 
develops during which the one-particle distribution relaxes towards the equilib- 
rium value with a relaxation time. This is characterized by a nonlocal kinetic 
equatioixin agreement with the virial correction to the approached equation 
of state Ej. We will focus on the first stage which is related to the formation 
of correlations and will find a measurable effect on the two-particle correlation 
function. The extracted lifetimes are shown to be too high if one ignores this 
transient time effects. This can be considered to belong to the initial state 
correlations. While inUO the initial state correlations are assumed to be of 
equilibrium type and discussed by density and temperature dependence we 
like to investigate here a nonequilibrium effect which consists in the fact that 
correlations need a specific time to be build up. 

The formation of correlations is connected with an increase of the kinetic 
energy or equivalently the build up of correlation energy. This is due to rear- 
rangement processes which let decay higher order correlation functions until 



2 



only the one - particle distribution function relaxes. 



2 Kinetic description 

The time dependence of the kinetic energy as a one-particle observable will 
be investigated within the kinetic theory. This can only be accomplished if 
we employ a kinetic equation which leads to total energy conservation. It is 
immediately obvious that the ordinary Boltzmann equation cannot be appro- 
priate for this purpose because the kinetic energy is in this case an invariant 
of the collision integral and constant in time. In ccmtLast, we have to consider 
non-Markovian kinetic equations of Levinson typellj'O, which account for t 
fnrrnatimi of two particle correlations and which conserves the total energy 

b 

t 

X j dt exp |-^| cos - t)\E^ [fjUah - fahfaK] , (4) 

to ^ 

where Ae — ^ h -rr ^^^^r-^ (p+i) denotes the energy difference between 

initial and final states. The retardation of distributions, fa{k,t), f'^{k — q,i) 
etc., is balanced by the quasiparticle damping Tg and / = 1 — /, the free particle 
dispersion E — /2m and the spin-isospin degeneracy s. The distribution 
functions are normalized to the density as J (^^fp/Cp) = n/s. 

The Boltzmann collision integral is obtained from equation if: (i) One 
neglects the time retardation in the distribution functions, i.e. the memory 
effects. This would lead to gradient contributions to the kinetic equation which 
can be shown to be Kespansible for the formation of .liigh energy tails in the 
distribution function EjO and virial corrections E2l'E3't3. This effect will be 
established on the second stage of relaxation, (ii) The finite initial time Iq 
is set equal to — oo corresponding to what is usually referred to as the limit 
of complete collisions. This energy broadening or off-shell behavior in is 
exclusively related to the spectral properties of the one-particle propagator and 
therefore determined by the relaxation of two-particle correlation. 

Since we are studying the very short time region after the initial dis- 
turbance we can separate the one-particle and two-particle relaxation. On 
this time scale the memory in the distribution functions can be neglected 
fait) = /a(0), exp I^T^I — 1 but ^'^^^ keep the spectral relaxation implicit 
in the off-shell cos- function of (||). 
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The eq. (^) can be integrated with respect to time and the resulting 
equation for f{t) represents the deviation of Wigner's distribution from its 
initial value, fa{t) = /a(0) + Sfa{t), and reads 

This formula shows how the two-particle and the single-particle concepts of the 
transient behavior meet in the kinetic equation. The right hand side describes 
how two particles correlate their motion to avoid the strong interaction regions. 
Since the process is very fast, the on-shell contribution to 6 fa, proportional to 
t/r, can be neglected in the assumed time domain and the Sf has the pure off- 
shell character as can be seen from the ofF-shcll factor A^^ (l — cos {j-IAe})- 
The off-shell character of mutual two-particle correlation is thus reflected in 
the single particle Wigner's distribution. 

The very fast formation of the off-shell contribution to Wignenls xListribu- 
tion has been found in numerical treatments of Green's functions E3'l3. Once 
formed, the off-shell contributions change in time with the characteristic time 
T, i.e., following the relaxation (on-shell) processes in the system. Accordingly, 
the formation of the off-shell contribution signals that the system has reached 
the stateJJre evolution of which can be described by the nonlocal Boltzmann 
equationll3, i.e., the transient time period has been accomplished. 

From Wigner's distribution (||) one can readily evaluate the increase of the 
kinetic energy 



ab 

This expression holds for general distributions fa- We choose a model of two 
initially counter - flowing streams of nuclear matter and plot in figure Q the 
time evolution of the kinetic energy and the correlation. 

This build up of correlations is found to be independent of the form of 
the initial distribution. If for example we choose a (equilibrium) Fermi dis- 
tribution as the initial distribution a build up of correlations will occur as 
well. This is due to the fact that the spatial correlations relate in momentum 
space to excitations, resulting in a distribution looking somewhat like a Fermi- 
distribution but wijLh a temperature higher than that of the initial uncorrelated 
Fermi-distributiontJ 
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Figure 1: The time evolution of the correlation and kinetic energy for two initially counter - 
flowing streams of nuclear matter with density and temperature ni = no/60 Ti = 0.5 MeV 
and n2 = no/20 T2 = 0.1 MeV moving with relative velocity of Ih/ fm which corresponds 
to a collision energy of 21 MeV/n. In this case the correlation time = 9/m/c. The 
equilibrated temperature is here (neglecting correlations) T = 11.6 MeV and a Gaui3 type 
of potential has been used. 



3 The formation of correlations 

We start now to calculate the observed increase of kinetic energy analyti- 
cally. As a first example we shall consider a Yukawa-type potential of the 
form Vyir) = 2^e~'"' where in nuclear physics applications gab is the cou- 
pling and K the effective range of potential given by the inverse mass of in- 
terchanging mesons. As a second example we shall use a Gaufi-type potential 
Voir) — Vbe"^''/''-' which has been used in nuclear physics application£HlElil 
with r/ = 0.57fm~^ and Vb = -453 MeV. 

The high temperature limit ol-th£ time dependence of the correlation en- 
ergy can be calculated analyticallyOEj. The low temperature value is of special 
interest, because it leads to a natural definition of the build up of correlations. 
One obtains 

d 1 V'^ ~ 

ut 2 cos I 

with standard notation of angular integrals and finds the time dependence 
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of the correlation energy E 



1 - - sin(a;) + ( 



X V T^T J \ 6 X 




3' ttT 
1 



sin(a;) 



(8) 



with X — 2^^j^t and the-epuihbriuin correlation energy E^°^^ for the Yukawa- 
or Gaufi- potential reads t2l 

f arctan ^ + ] Yuk. 




^erf(^) GauP 

(9) 

where e/ = -^^^ and hi = The best choice for cut-off we found Ec « 

£/. We see that the correlation energy is built up and oscillates around the 
equilibrium value damped with t^^ in time. We now define the build up time 
Tc as the time where the correlation energy has reached its first maximum. 
This time is given by Tc « 

We assert here that the result (||) is valid for any binary interaction. We 
could use Born or T-matrix approximation and the same time dependence but 
different E]"^ would result. This correlation time limits the ■validity of quasi- 
particle picture which is established at times greater than Tc Ej. Incidentally, 
in the early 1950s the criterion h/ksT < r mas supposed to limit the validity 
of the Landau Fermi-liquid theory for metalsE3. Later it was shown by Landau 
that this criterion is irrelevant and he proposed the correct criterion t > h/ep. 

In order to illustrate the weak temperature dependence of Tc we plot in fig- 
ure H (thick lines) results from the solution of the Kadanoff and Baym equations 
for a fixed chemical potential of 37.1 MeV and for three different temperatures. 
The figure shows the increase of the kinetic energy (equivalent to the decrease 
of correlation energy) with time. The KB results are compared with those 
from approximation (||). One sees that the agreement is good initially while 
correlations are built up. At low temperatures the oscillations discussed above 
are obvious in the approximate results while the KB calculations only show a 
slight overshoot at the lowest temperatures. We believe that the discrepancy 
is due to the damping that is neglected in (^. 

We can give a Pade formula for the increase of kinetic energy which inter- 
polates the results of K/B equation within the temperature range up to Fermi 
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Figure 2: The formation of correlation plotted as an increase of the kinetic energy with time 
for temperatures 1,10,40 MeV. The chemical potential is fixed to 37.1 MeV which corre- 
sponds to densities 0.16, 0.18, 0.35fm~^ . The thick lines show results from KB calculations 
while the thin lines are approximate values via formula feh. The equilibrium correlation 
energy was chosen to be equal to the KB result. The oscillations are overestimated by the 
apprnptipiate formula. For T = 40 MeV we plotted also the high temperature approximate 
value as thin solid line. The built up of correlations is too fast. 



energy Ck/b and the analytical result (||) signed ({t) 

_ Ecorrjt) _ 1.5225 (1 - cos(0. 7337a:) + sin(0.6568x)) 

at) = 1 - (10) 
X 

with X — Aept/fi. The Ck/b is a smooth curve through C,{t). As we pointed 
out above, this result is universal for any binary interaction approximation. 
Therefore we can adopt this scaling as a universal scaling of the one-particle 
Wigner distribution. Provided we have solved any local Markovian kinetic 
equation like BUU, we can simply multiply the obtained Wigner function by 
the time dependent factor ([l0| ) in order to incorporate the transient time effect 
approximately. Since we fixed this factor to the observable kinetic energy, we 
expect to have correctly described the energy variables. 



7 



4 Application to interferometry 

Now we employ the found scaling of short time effects to the formula and 
adopt a space-time source parameterization of 

.g(p, r, t) = poe(r)e(i?, - r)e(i)pe-^e"*/^CW (H) 

with the radius Rs and lifetime r of the source and the scaling ({t) of (p^. 
Without center of mass momentum P the correlation function becomes-.inde- 
pendent of lifetime and collapse to the standard interferrometry resultl^J 

Considering the momentum dependence and assuming initially free particles 
\(p{p,r)\'^ oc 1 ± cos(rg) we obtain a separable representation 

R{P, q) + l = l± ^{q.RsMC, q.p, r) (13) 

with 

= ((-3 + 3x^ - 2x^ + Ax^) cos{2x) + (-6x - &x^ - 4a;'*) sin(2a;)) 
2x° 

(14) 

and 

oo 

/ didt'e-(*+*')/"C(OC(^') cos{q.p{t' - t)/2m) 
n(C, (Z.p/2m, r) = — . (15) 

While the first factor 5* gives the dependence of the correlation function 
on the source size, the second factor contains the effect of finite lifetimes. 
The result without transient time effects is that with increasing lifetime the 
correlation functions are diminishedLi. One gets for C = 1 just 

"°(^''^-^'")= l + (Jr/2m)^ - 
We can write the analytical result for (^ with C,{t) from ( p^ 

n(a, c, r) = I 4 ( — ^— + arctan((a — c) t) + arctan((a + c) t) 
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Eigjure 3: The ratio of correlation funetion 11 of ( 
(ll™ without short time dynamics versus hfetime of 



) with short time dynamics to Ho of 
e source. 



log(l + (a - c)^ t2) + log(l + (a + c)^ r^) 



2 



/ (16 (ar — arctan(aT)) ) 



(17) 



with c = Atp/h and a = q.p/2m. Figure y shows the influence of the transient 
time effect on the correlation function. We have plotted conveniently the ratio 
of reduced correlation functions 11 of ( |l5|) with short time dynamics to ITq of 
( |l6| ) without short time dynamics. 

We recognize that there is a suppression for smaller lifetimes. Since the 
effect of lifetime r itself lowers the correlation function we can state that the 
extracted lifetimes from experiments arc too high if short time or transient time 
effects are neglected. For higher lifetimes this effect is opposite. The extracted 
lifetimes from experiment should be too low. Interesting is the density depen- 
dence. While first with decreasing density the lowering/enhancement shifts to 
higher lifetimes, there is a critical density where we have only suppression. It 
is found to be around 0.0036fm~'^. 

Due to the found universal scaling this effect is important in any recent 
simulation of BUU- type since the transient time effect of formation of cor- 
relations and therefore increase of kinetic energy has not been considered so 
far. For moleculardynamical methods this is in principle included but buried 
in artificial initial correlations due to the nuinerical set up. For a comparison 
between the kinetic and MD simulation seella. 
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5 Summary 



The gradient approximation of the kinetic equation in second order Born ap- 
proximation is investigated. A finite duration approximation of the non Marko- 
vian colhsion integral is proposed which follows from time dependent Fermi's 
Golden Rule and which is in good agreement with the numerical solution of 
complete Kadanoff and Baym equation. 

The build up time of correlations is investigated and it is found that the 
low temperature value is universal for any approximation at the binary collision 
level. It is shown that the formation time of correlations is nearly determined 
by the ratio of h to the transfer energy which can be considered as an analogue 
to the uncertainty principle. 

A universal scaling is formulated which allows to map the Wigner function 
of Markovian simulation to a result including the transient time effect. As a 
consequence the influence to interferrometry methods is discussed. It is found 
that the extracted lifetimes from experiment should be too large for lower 
lifetimes and too small for higher values depending on the freeze out-density 
if transient time effects are neglected. 

Acknowledgments 

The authors like to thank P. Lipavsky and V. Spicka for interesting discussions 
and N. Kwong for many useful hints. 

References 

1. G. F. Bertsch, P. Danielewicz, and M. Herrmann. Phys. Rev. C, 49:442, 
1994. 

2. D. H. Boal, C. K. Gelbke, and B. K. Jennings. Rev. Mod. Phys., 62:553, 
1990. 

3. R. Hanbury-Brown and R. Q. Twiss. Nature (London), 178:1046, 1956. 

4. J. Pochodzalla and et. al. Phys. Rev. C, 35:1695, 1987. 

5. S. E. Koonin. Phys. Lett. B, 70:43, 1977. 

6. W. G. Gong and et. al. Phys. Rev. C, 43:1804, 1991. 

7. W. Bauer. Nucl. Phys. A, 471:604, 1987. 

8. S. Pratt and M. B. Tsang. Phys. Rev. C, 36:2390, 1987. 

9. W. G. Gong, W. Bauer, C. K. Gelbke, and S. Pratt. Phys. Rev. C, 
43:781, 1991. 

10. L. P. Kadanoff and G. Baym. Quantum Statistical Mechanics. Benjamin, 
New York, 1962. 

11. H.S. Kohler. Phys. Rev. E, 53:3145, 1996. 



10 



12. N. N. Bogoljubov. J. Phys. (USSR), 10:256, 1946. transl. in Studies 
in StaUstical Mechanics, Vol. 1, editors D. de Boer and G. E. Uhlenbeck 
(North-Holland, Amsterdam 1962). 

13. M. Bonitz and et. al. J. Phys.: Condens. Matter, 8:6057, 1996. 

14. P. Lipavsky, F. S. Khan, A. Kalvova, and J. W. Wilkins. Phys. Rev. B, 
43(8):6650, 1991. 

15. K. Morawetz, V. Spicka, and P. Lipavsky. Formation of correlations in 

plasma at short time scales. Phys. Lett. A in press. 

16. V. Spicka, P. Lipavsky, and K. Morawetz. Phys. Lett. A, 240:160, 1998. 

17. T. Aim. PhD thesis, Universitat Rostock, 1992. report GSI-93- 
07, Darmstadt. 

18. T. Aim, G. Ropke, and M. Schmidt. Phys. Lett. B, 301:170, 1993. 

19. L B. Levinson. Fiz. Tverd. Tela Leningrad, 6:2113, 1965. 

20. L B. Levinson. Zh. Eksp. Tear. Fiz., 57(2):660, 1969. [Sov. Phys.- 
JETP 30, 362 (1970)]. 

21. K. Morawetz. Phys. Lett. A, 199:241, 1995. 

22. A. P. Jauho and J. W. Wilkins. Phys. Rev. B, 29(4):1919, 1984. 

23. K. Morawetz, R. Walke, and G. Ropke. Phys. Lett. A, 190:96, 1994. 

24. H. Haug and A. P. Jauho. Springer, Berlin Heidelberg, 1996. 

25. K. Morawetz and G. R5pke. Phys. Rev. E, 51(5):4246, 1995. 

26. V. Spicka and P. Lipavsky. Phys. Rev. B, 52(20):14615, 1995. 

27. V. Spicka, P. Lipavsky, and K. Morawetz. Phys. Rev. B, 55(8):5084, 
1997. 

28. V. Spicka, P. Lipavsky, and K. Morawetz. Phys. Rev. B, 55(8):5095, 
1997. 

29. P. Danielewicz. Ann. Phys. (NY), 152:305, 1984. 

30. H. S. K5hler. Phys. Rev. C, 51:3232, 1995. 

31. H.S. K5hler. Phys. Rev. C, 51:3232, 1995. 

32. P. Danielewicz. Ann. Phys. (NY), 152:239, 1984. 

33. K. Morawetz and H. S. Koehler. Phys. Rev. C, 1997. sub. 

34. H. Smith and H. Hojgaard-Jensen. Clarendon Press, Oxford, 1989. 

35. R. Peierls. Quantum Theory of Solids. Oxford University Press, London, 
1955. 

36. B. K. Jennings, D. H. Boal, and J. C. Shillcock. Phys. Rev. C, 33:1303, 
1986. 



11 



